Structure and stability of small H clusters on graphene 
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The structure and stability of small hydrogen clusters adsorbed on graphene is studied by means 
of Density Functional Theory (DFT) calculations. Clusters containing up to six H atoms are inves- 
tigated systematically - the clusters having either all H atoms on one side of the graphene sheet 
( CIS-clusters) or having the H atoms on both sides in an alternating manner (trans-cluster). The 
most stable cis-clusters found have H atoms in ortho- and para-positions with respect to each other 
(two H's on neighboring or diagonally opposite carbon positions within one carbon hexagon) while 
the most stable trans-clusters found have H atoms in ortho-trans-positions with respect to each 
other (two H's on neighboring carbon positions, but on opposite sides of the graphene). Very stable 
trans-clusters with 13-22 H atoms were identified by optimizing the number of H atoms in ortho- 
trans-positions and thereby the number of closed, H-covered carbon hexagons. For the cis-clusters, 
the associative H2 desorption was investigated. Generally, the desorption with the lowest activation 
energy proceeds via para-cis-dimer states, i.e. involving somewhere in the H clusters two H atoms 
that are positioned on opposite sites within one carbon hexagon. H2 desorption from clusters lacking 
such H pairs is calculated to occur via hydrogen diffusion causing the formation of para-cis-dimer 
states. Studying the diffusion events showed a strong dependence of the diffusion energy barri- 
ers on the reaction energies and a general odd-even dependence on the number of H atoms in the 
cis-clusters. 
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I. INTRODUCTION 



The interaction of H atoms and molecules with car- 
bon based materials like graphite, single wall carbon nan- 
otubes (SWCNT) and graphene has attracted consider- 
able interest during the last two decades, since these sys- 
tems are of interest in fields as diverse as astrochemistry, 
hydrogen storage and nanoelectronics. 
The understanding of H chemisorption on graphitic ma- 
terials and formation of H2 molecules at low H atom den- 
sities is highly relevant for interstellar chemistry!,. H2 
formation in interstellar dust and molecular clouds is ex- 
pected to occur via recombination of H atoms adsorbed 
on interstellar dust grain surfaces. Since carbonaceous 
grains are abundant in the interstellar medium, particu- 
lar attention is focused on H chemisorption on graphite 
surfaces^. 

Most of the recent efforts regarding H adsorption on 
carbon-based materials have been stimulated by the pos- 
sibility of using carbon nanostructures as a hydrogen 
storage mediunt^. Earlier reports on H storage of ~7.4 wt 
% in nanostructured graphite^ and 14 wt % in SWCNTs^ 
have not been confirmed by new experiments that rather 
indicate that the H storage capacity at room tempera- 
ture is 3.8 wt % in graphite nanofiberSi^ and ~5.1 % in 
SWCNTs^. The precise mechanism of H adsorption in 
these structures is unknown since these values are well 
above the estimated limit of 1% for H2 physisorption at 
room temperature-. For graphene a theoretical maxi- 
mum storage capacity of 7.7 wt % is reported-. 
Recent fabrication of graphene has generated an enor- 
mous number of new studies focused on the prospects 
of using graphene as a key material in post silicon 



electronics^"—. The high electron mobility can be ex- 
ploited for construction of ballistic transistors operating 
at frequencies unreachable with current semiconductor 
materialaiS. Yet, for applications in nanoelectronics effi- 
cient methods and tools for tailoring the electronic prop- 
erties of graphene are required. One of the most promis- 
ing routes towards the opening of a band gap in graphene 
is based on its hydrogenation^ii^^J^. Hence, an atomic 
scale description of the binding mechanism and stability 
of small hydrogen structures formed on graphene will en- 
hance efforts aimed at full control over graphene band 
gap engineering with deposited H atoms. 
The chemisorption of H monomers on graphite has 
been studied extensively both experimentally and 
theoreticallyi^"-^^. H dimers on graphite or graphene 
have also been thoroughly studied in recent yearsi^r^. 
According to experiments, bigger structures will form 
when graphite or graphene are exposed to higher H 
doses^^i^>^>^>^. These structures, highly relevant for 
possible technological applications, are much less inves- 
tigated than dimers. Only recently, Casolo et al.^*, Ro- 
man et alM, Ferro et al^ and Khazaei et al^ used 
Density Functional Theory (DFT) to calculate the en- 
ergetics of hydrogen clusters with three and four atoms, 
as well as one structure with six atoms on a graphene 
sheet. Cuppen et al^ studied the formation of hydro- 
gen clusters using Kinetic Monte Carlo methods. Luntz 
and co-workerSi22, studied desorption of D2 molecules 
from D clusters adsorbed on graphite, combining ex- 
perimental and theoretical methods. A comprehensive 
list of tetramer structures and the recombination path- 
ways was produced based on measured desorption den- 
sity and DFT calculations'^'^. In the present paper we 
apply DFT to systematically study small hydrogen clus- 
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ters composed of three to six H atoms on graphene. Their 
structure and stability against H diffusion and H2 recom- 
bination are determined at the atomic scale. Two types 
of clusters are considered, czs-clusters having all H atoms 
on the same side of the graphene sheet, and trans-clusters 
having the H atoms on both sides of the graphene sheet 
in an alternating manner. For the trans-clusters a few 
larger clusters were further considered. The computa- 
tional method is described in Section |lTl The calculated 
H structures and their energetics are presented in Section 
mil The discussion of the results, including observed 
trends in the H binding with the size of H clusters, is 
given in Section ITVl The results are summarized in Sec- 
tion |Vl 
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TABLE I. Binding energies (in eV) of nH (1< n < 6) cis- 
clusters adsorbed on graphene, as a function of the simula- 
tion cell size and the number of k-points used to sample the 
Brillouin zone. Results are given for the monomer, the ortho- 
dimer ( "O" ) and the most stable structures from Table IIIII 



II. COMPUTATIONAL DETAILS 

The DFT calculations were performed with the plane 
wave based DACAPO program package^i^, applying 
ultra-soft pseudopotentials'^^''^'' to describe electron-ion 
interactions, and the Perdew Wang functional (PW91) 
for the electronic exchange correlation effects. The elec- 
tron wave functions and augmented electron density were 
expanded in plane waves with cutoff energies of 25 Ry 
and 140 Ry, respectively. The H cluster configurations 
were calculated by modelling the graphite surface with 
rhombohedral, periodically repeated slabs consisting of 
one graphene sheet having a 6x6 surface unit cell with 
72 carbon atoms, and separated by 15 A of vacuum. The 
Chadi-Cohen scheme^ with six special points was used 
for sampling of the surface Brillouin zone. Binding en- 
ergies of the nH clusters are reported per cluster using 
the clean graphene sheet plus n separate H atoms as the 
reference system, i.e.: 

E^, = {E[gra) + n£;(H)] - £:(nH/gra), (1) 

where EinWj gra)^ E{gra) and i?(H) are total energies 
of the TiH clusters at graphene, the clean graphene sheet 
and the free H atom, respectively. The binding energies 
are well converged with respect to the number of k-points 
and the surface cell size as evidenced by Table |T] present- 
ing the binding energies of the most favorable nH cis- 
clusters (n < 6) on graphene that we have encountered 
in the present investigation (i.e. the monomer, the ortho- 
dimer "O" , and trimers up to hexamers to be defined 
below). The H binding energy in gas-phase H2 molecule 
calculated using computational method described above 
is 2.29 eV/H. 

In figures [T][7] schematic illustrations of the H islands 
omitting the ionic relaxation patterns are show. How- 
ever, all structures were fully relaxed using the Broyden- 
Fletcher-Goldfarb-Shanno algorithmSE. The activation 
energies for diffusion of H atoms and for the H2 recombi- 
nation were calculated applying the nudged elastic band 
method^ using at least seven configurations to model re- 
action paths. For sequences of reaction steps, we use the 
term net barrier to denote the energy difference between 



the initial state and the highest potential energy point 
along the path. This highest energy point will be a tran- 
sition state of some n'th reaction step along the path and 
the net barrier becomes the sum of the reaction energy 
for moving from the first initial state to the n'th initial 
state plus the local barrier in the n'th reaction step. 



III. RESULTS 

In the present paper we considered two types of ad- 
sorbed nH clusters: (A) cis-clusters with all H atoms ad- 
sorbed on the same side of the graphene sheet, and (B) 
trans-clusters where adsorbates bind at both sides of the 
graphene in an alternating manner. The corresponding 
results are presented in subsections IIII A l and IIII 5] 



A. Cis-clusters: H adsorbed only on one side of the 

sheet 

Since the activation energy of H diffusion through 
graphene is higher than 4 eV*^ in the majority of ex- 
periments related to the H adsorption on graphite or 
graphene, the adsorbates will exclusively bind on the side 
of the graphene sheet exposed to the source of H atoms. 
Hence, we first consider this class of H configurations on 
graphene. 

1. Hydrogen monomer and short dimers 

The structure and binding energies of the H 
monomer and dimers have been studied by several 
groupsi^"— i^iiS^iiSiii^. Our recent publication^^ contains 
well tested results for the H monomer and a comprehen- 
sive list of hydrogen dimer structures on graphene. The 
study includes H binding energies and activation ener- 
gies of the most important kinetic processes. Since the 
results for bigger clusters presented in this paper will be 
compared to those of monomers and short dimers, for 
the sake of convenience we include in Fig. [T] and Table [Tl] 



3 



monomer dimers 



config. binding 


config 


binding 


energy 




energy 


monomer 0.77 





2.73 




P 


2.68 




M 


1.58 



TABLE II. Binding energies (in eV) of tlie monomer and sfiort 
dimers at graphene^^. 

relevant results of Ref<^, calculated with the same com- 
putational set-up as used in the present study. 




FIG. 1. (a) Three H short-dimer configurations on graphene. 
H and C atoms are represented with small black and gray 
spheres, respectively. Energies are given relative to the energy 
of the O-dimer structure; (b) the path and potential energy 
diagram for H2 recombination from the O-dimer configura- 
tion. The transformation from the O-dimer to the P-dimer 
has a net barrier of 1.65 eV (the sum of 1.15 and 0.50 eV). 
All energies are in eV. 



2. Hydrogen trimers 

The trimer structures investigated in the present work 
are depicted in Figure [DJa), together with the total ener- 
gies calculated relative to the energy of the Ha^ structure, 
identified as the most stable among the trimers consid- 
ered. The binding energies of the H trimers are given 
in Table IIIII It is instructive to interpret the three 
most stable trimers (i.e. Hs^, H^^^ and Ha^'^) as being 
built from an ortho- or para-dimer (O- or P-dimer) with 
an extra atom added in an O- or a P-position. Adding 
the extra atom a in meta-position (M-position) results 
in less stable structures (i.e. Ha^^ and Ha^). Thus the 
presence of an embedded M-dimer, i.e. involving some- 
where in the H clusters two H atoms in M-position with 
respect to each other, appears to be strongly disfavored. 
The least stable trimers considered are triangular struc- 
tures with the atoms in M- or extended dimer-positions 
(i.e. Ha HaV" and Ha^i"). 

The Ha^ and Ha"i structures were also presented by 
Casolo et al^, with the binding energies in very good 




H2(gas) + H2(gas) + 

monomercagraphene monomer@graphene 



FIG. 2. (a) H trimer configurations on graphene. Energies 
are given relative to the energy of the structure; (b) the 
paths and potential energy profiles for H2 recombination from 
the H3I and Hs^' configurations. All energies are in eV. 



agreement with the values reported here. The Ha 
trimer was suggested by Khazaei et al^ as a model for 
experimentally observed star like STM patterns, but ap- 
pears as one of the least favorable trimer structures ac- 
cording to our work. Since we further calculate an energy 
barrier for H diffusion to neighbouring C sites of less than 
0.65 eV for the Ha^^i configuration, it is unlikely to be 
stable at temperatures up to 600 K— . A very compre- 
hensive list of trimer structures is reported by Roman 
at al^. The most favorable structures ( Ha^, Ha^^ and 
Ha™) observed in the present work are in full agreement 
with the results in Ref.^^. 

Considering the associative desorption of H2 from the 
trimer structures we find that it occurs via an embed- 
ded P-dimer. The calculated pathways are sketched in 
Fig. m^b) for the two most favorable structures (i.e. Ha^ 
and Ha^i). For configuration Ha^ that already contains 
an embedded P-dimer, the barrier is calculated to 1.38 
eV, which is close to the one found for the isolated P- 
dimer (1.40 eV3^ - see Fig. ^h)). At the transition 
state along this pathway the H2 molecule is nearly par- 
allel to the surface with the H-H bond-length of ~1.2 A. 
Very similar geometries of transition states are accounted 
for other reaction paths shown in Figs. [3] and [5j The 
H3II trimer does not contain an embedded P-dimer and 
the associative desorption from the embedded O-dimer 
is calculated to be associated with a barrier of 1.61 eV. 
An alternative scenario includes several H diffusion steps 
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TABLE III. Binding energies (in eV) of nH cis-clusters at graphene. The cluster configurations are depicted in Figs. [2]-[5] 



transforming the H3 trimer via the H3 trimer to the 
trimer for which the associative desorption can occur. 
Along this path, the highest activation energy becomes 
that of the transition from the H3'' to the H3^ trimer, 
which is 1.49 eV, see Fig. [^b). 



3. Hydrogen tetramers 

The most favorable tetramer structures found are il- 
lustrated in Figure Ela) and their binding energies are 
given in Table IIIII There are many low energy tetramers 
and a common feature for the very most stable of them 
(i.e. H4'~^'') seems to be that they are composed of one 
of the most stable trimers with an extra atom attached 
in an O- or a P-position. It is interesting to note that a 
triangular arrangement of the atoms seems to be particu- 
larly unfavorable. Comparing H4''^' and £[4"^'' one might 
guess based on the previous discussion that H4-'^'' would 
be the more stable, since it contains exclusively atoms 
in P-positions, whereas in H4-^^ one atom is moved to a 
M-position. The same arguments apply when comparing 
H4^''^ and H4''^. These triangular structures are unfa- 
vorable due to three H atoms adsorbed on the same C 
sublattice. Such imbalance in the occupation of two car- 
bon sublattices creates additional unpaired electrons in 
the graphene lattice, increasesing the total energy of the 
system. Note that for the trimers considered triangular 
structures were also very unfavorable. The pathways for 
H2 recombination from the four most favorable configura- 
tions are depicted in Figure |3fb). As for the trimers, the 
recombination occurs from an embedded P-dimer, either 
directly [for the H4'\ H4''' and H4^^ tetramers], or via H 
diffusion from the H4' over the H4'^''' to the H4''' struc- 
ture, followed by H2 associative desorption, as shown in 
Figure Efb) . The net barrier for H2 formation from the 
H4^ tetramer becomes 1.81 eV (the sum of 1.44 eV and 
0.37 eV), which demonstrates a high stability of this con- 



figuration. Other possible scenarios for H2 recombination 
from tetramers are discussed in Refi^. 



4-. Hydrogen pentamers 

The most stable pentamer structures in Figure |4] are 
formed by attachment of an additional H atom in an O- 
or a P-position to one of the preferential tetramer struc- 
tures. The binding energies of the configurations con- 
sidered are provided in Table IIIII The binding energies 
of the seven most stable structures are within an energy 
window of 0.2 eV. Although the other studied configura- 
tions are less favorable, the energy difference between the 
least (the Hs-^'^) and the overall most stable structure 
(the II5') is smaller than 0.6 eV. 

The activation energies for associative desorption are not 
calculated, since we expect that the H2 formation occurs 
along pathways that are qualitatively similar to those of 
the tetramers, i.e. via direct associative recombination 
from embedded P-dimers or several steps of H diffusion 
between adjacent C sites, prior to the recombination from 
a P-dimer state. 



5. Hydrogen hexamers 

Hexamers are the largest H clusters considered sys- 
tematically in the present study. In Fig. [SJa) the most 
stable structures identified are shown. The most favor- 
able structure, Hg^ has a total binding energy of 9.62 
eV. Following the general trend outlined, the most stable 
hexamers contain motifs from stable smaller clusters and 
have all atoms in O- or P-positions. It is interesting to 
note that the Hg^ and the Hg^' structures have almost 
identical stability despite being quite different in geom- 
etry; the Hg^ is of low symmetry and has H atoms in 
0-positions while the Hg^' is of high symmetry and has 
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H2(gas) + P 

FIG. 3. (a) H tetramer configurations on grapliene. Potential 
energies are given relative to the energy of the H4' structure, 
(b) The pathways for H2 recombination from the H4', H4^\ 
and H4^^ configurations. The OPO is one of the long- 
dimer configurations from Ref»Si,. All energies are in eV. 



H atoms in P-positions. Generally, very many hexamer 
structures are found that differ only by a few tenths of 
an eV. For most of the structures presented in Fig. ^a.) 
the mechanism for the H2 recombination and correspond- 
ing energy barriers are expected to be very similar to 
those determined for the dimers, trimers and tetramers. 
Thus, we did not perform the corresponding calculations. 
The reaction paths were however carefully calculated for 
the He^^ and Hg^^^^ structures, since both configurations 
possess ideal hexagonal symmetry. The breaking of the 
high symmetry of these two configurations might lead to 
activation energies for H2 recombination that are higher 
than the typical values found for other H structures. The 
results for the H2 formation from the Hg^^ and Hg^^^^ 




FIG. 4. H pentamer configurations on graphene. Potential 
energies (in eV) are given relative to the energy of the H5' 
structure. 

hexamers are presented in Fig. [SJb). Indeed, the activa- 
tion energies of 1.72 eV and 1.88 eV, calculated for Hg^^ 
and Hg^^^^, respectively, indicate high kinetic stability of 
these structures with respect to molecular hydrogen for- 
mation. Ferro et al.^- have suggested the Hg^^ hexamer 
structure as a possible candidate for the experimentally 
observed very stable starlike STM patterns^l. 



B. Trans-clusters: H adsorbed on both sides of the 

sheet 

A common feature of the dimers in ReL^^ and the cis- 
clusters presented in Figures [5] to [5] is that all H atoms 
are adsorbed on the same side of the graphene sheet. 
Given a very high energy barrier for H diffusion through 
the graphene layer this scenario is most likely to occur in 
defect-free samples deposited at crystalline surfaces. Yet, 
in free-standing graphene, graphene samples with defects 
or in small graphene patches, diffusion barriers at the 
defects or edges could be significantly lower than at the 
perfect sheet, opening routes for H adsorption on both 
sides of the layer. It turns out that such clusters - trans- 
clusters - are much more stable than those considered 
in the previous section. In the following, we present the 
results of our systematic investigation of trans-clusters 
with 2-6 H atoms, showing, however, only the very most 
stable structures identified. We further selectively chose 
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FIG. 5. (a) H hexamer configurations on graphene. Potential 
energies are given relative to the energy of Ho' structure, (b) 
The reaction paths for H2 recombination from the He'' and 
He 



configurations. All energies are in eV. 



to study a few larger trans-clusters that will be presented 
and discussed. 

The 0-trans-dimer formed by adsorption of two H atoms 
in the configuration depicted in Fig. is very stable 
with a total binding energy of 3.30 eV. This is by 0.54 
eV higher than the binding in the 0-dimer with the H 
atoms on the same side of the graphene sheet, i.e. the O- 
cis-dimer. The high stability of the O-trans-dimer con- 
figuration has already been reported by Roman et al.^'^. 
The high stability of the O-trans-dimer appears to be re- 
lated to the immediate proximity of the C atoms, since 
expanding the trans-dimer forming the P-trans- and M- 
trans-dimers of Fig. and Fig. [6j; leads to binding en- 
ergies that are smaller by 0.79 and 1.82 eV, respectively, 
than that of the O-trans-dimer. In view of this finding, 
the search for the structures of the most stable larger H 
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FIG. 6. H structures with atoms adsorbed on both sides of 
the graphene sheet: (a)-(c) trans-dimers; (d) trans-trimer; 
(e)-(g) trans-tetramers; (h)-(i) trans-pentamers; (j)-(k) trans- 
hexamers. Adatoms from opposite sides of the sheet are col- 
ored differently. Binding energies per H atom are given in 
eV. 



clusters was limited to structures in which all neighbor- 
ing H atoms were configured as in the O-trans-dimer. 
The Hs^^-trans-trimer structure shown in Fig. [6ji is the 
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TABLE IV. Binding energies of hydrogen trans-clusters ad- 
sorbed on graphene. 

most favorable trans-trimer structure found in our study. 
The binding energy of 5.11 eV, is by 0.89 higher than in 
the Ha^ cis-trimer. 

Two almost equally stable trans-tetramers (H4^-trans 
and H4^-trans) were found as shown in Fig. |6j3-f with 
binding energies of 7.56 and 7.36 eV. The triangular 
shaped trans-tetramer in Fig. |6^ is slightly less stable 
with a binding energy of 6.83 eV, which is 0.73 eV smaller 
than Fig. IBJ3. This indicates that also for trans-clusters 
a triangular arrangement of the H atoms is disfavored. 
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and can in the same way as for the triangular shaped 
cis-clusters be related to an imbalance in the number of 
H atoms adsorbed on the two C sublattices. The most 
stable configuration, Fig.|B}3, is as much as 1.26 eV more 
stable than the H4^-cis-tetramer in Fig. [3l 
Two H trans-pentamers (Hs^-trans and Hs-^^^-trans) were 
investigated as shown in Figl^h and [SJ. The calculated 
total H binding energies were 9.44 and 9.20 eV, respec- 
tively. The structure in Fig.|6)i is by 1.71 eV more stable 
than the Hs^-cis pentamer in Fig. U) 
The trans-hexamer in Fig. |B] (Hg^^^^-trans) is a particu- 
larly stable structure with a total binding energy of 12.47 
eV which is 0.66 eV larger than that of the configuration 
in Fig. [6ji (He^-trans), and as much as 2.79 higher than 
the binding of the most stable cis-structure, Hg^, in Fig. 

m 

The graphane, fully hydrogenated graphene with hy- 
drogen atoms adsorbed at carbon atoms on both sides 
of the sheet in an alternating manner^^ii^, can be consid- 
ered as the infinite trans-cluster. The H binding energy in 
graphane is 2.49 eV per atom, which is considerably more 
than the 2.08 eV calculated for the Hg^^^^-trans hexamer. 
Thus, with an increase in the size of trans-clusters we ex- 
pect higher H binding energies, which should approach 
the value calculated for graphane. To examine this ex- 
pected trend we considered several bigger trans-clusters. 
The larger trans-clusters considered are shown in Fig. 
[71 They have all been constructed as truncated pieces 
of graphane embedded in graphene. The most stable of 
the selected clusters are exclusively composed of closed 
hydrogenated carbon hexagons, already identified as the 
most stable small trans-clusters (Hg^^^^-trans). To fur- 
ther study the stability of these structures we included 
in our investigation configurations with missing H atoms 
on carbon hexagons at the cluster edges, which results in 
the occurrence of embedded M-dimers. Apart from the 
clusters having various sizes the primary difference is the 
configuration of the edge atoms. The H binding energies 
quoted in the figure do indeed increase with cluster size. 
This will be analyzed further in the Discussion Section. 



C. Hydrogen induced magnetism 
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In a bipartite lattice, any imbalance in number of 
sites belonging to each of two sublattices leads to a 
magnetic ground state^iii^ii^. One of the possible 
mechanisms to induce such an imbalance on graphene 
is adsorption of H atoms as has been discussed in the 
recent Iiterature22i2iii2ii^. According to our calculations, 
upon adsorption of a single H atom on graphene, small 
magnetic moments are observed at several C atoms in 
its vicinity. The calculated magnetic moments at the 
individual atoms are smaller than 0.1 /ib- Concerning 
H dimers, non-zero spin density is observed only for 
structures with both hydrogens adsorbed on the C 
atoms from the same sublattice. According to the DFT 
calculations, these structures are unfavorable, and there- 



FIG. 7. Bigger H structures with atoms adsorbed on both 
sides of the graphene sheet. Adatoms from opposite sides of 
the sheet are colored differently. Binding energies per H atom 
are given in eV. 



fore not expected to form in experiments. The energy 
gain due to spin-polarization does not affect the relative 
stability of the most favorable dimer structures obtained 
from non-spin-polarized calculations. This has already 
been demonstrated in our previous publicatior*2£. A 
similar effect of the spin-polarization is found in the 
present work for the energetics of the investigated H 
cis configurations with three, four, five and six atoms. 



a) 


c, 


c 


Id) 0)0 






12 3 4 5 6 7 



number of H atoms 

FIG. 8. The variation of the H binding energy with the size 
of cis- (circles) and trans-clusters (squares). The values cor- 
responding to the most stable cis-structures (O-dimer, H3'- 
trimer, H4'-tetramer, Hs'-pentamer and He'-hexamer) are 
shown as big circles. Results for less stable cis-clusters from 
Tables Ull and Hill are given as small open circles. 



However, several favorable trans-clusters investigated 
here (H4"^-trans, Hi2^^-trans, and H22^-trans) carry a 
total magnetic moment of 2 /ib- We note that for these 
clusters, there is an imbalance between the number of H 
atoms on the two sides of the graphene sheet and that 
such graphane islands embedded in the graphene lend 
themselves as possible building blocks for graphene- 
based magnetic materials. 



IV. DISCUSSION 

The H adsorption on graphene in compact cis-cluster 
structures, described in Sec. IIIIl is thermodynamically 
preferential compared to the adsorption of isolated H 
atoms. The binding energy of the H monomers is 0.77 
eV. Already for the smallest clusters, such as 0-dimers 
and trimers the binding per H atom increases 

to 1.37 and 1.39 eV, respectively. The H binding in 
tetramers is 1.57 eV, and slightly smaller in Hr,^ 
pentamers, where we calculated binding energy of 1.53 
eV per H atom. The strongest binding is evaluated for 
Hg^ hexamers, with the value of 1.60 eV per H atom. 
The evolution of the H binding energy with the size of 
the clusters is depicted by circles in Fig. [H] 
If we consider the binding configurations in compact 
trans-cluster structures with H atoms adsorbed on both 
sides of the graphene sheet, Fig.[6l the trend in calculated 
values is similar: 1.65 eV (dimer), 1.7 eV (trimer), 1.89 
eV (tetramer and pentamer), and 2.08 eV (hexamer). 
These binding energies are depicted by squares in Fig. [5] 
and are seen to be shifted to significantly higher values 
than in the clusters with H atoms adsorbed only on one 
side of the graphene sheet. Although our geometry for 
the 0-trans dimer is slightly different from the one found 
by Boukhvalov et al.^, their explanation of particularly 
favorable H adsorption on graphene observed in 0-trans 
dimers is fully applicable to the clusters considered in 
this study. In Fig. [H^-b side views of the 0-cis- and 



FIG. 9. The side views of (a) O-cis-dimer; (b) O-trans- 
dimer; (c) He^'^'-cis-hexamer; (d) He^'^'-trans-hexamer. The 
C atoms decorated with H are denoted as Co and their near- 
est neighbours without adsorbed H atoms as Ci. Adatoms 
from opposite sides of the sheet are colored differently. 



0-trans-dimers are shown and the carbon atoms are 
labeled. Using these labels, the H-C-C and C-C-C angles 
in the structures are given in Table |Vl From the Table 
it is seen that the angles calculated for trans-clusters are 
closer to the ideal tetrahedron value of 109.5° than those 
resulting for the cis-structures. Thus, in trans-clusters 
the H binding configurations allow creation of nearly 
perfect tetrahedral surroundings of C atoms directly 
involved in the interaction with adsorbates which results 
in an additional gain in the H chemisorption energy. An 
effect similar to that identified for the 0-dimers is found 
when comparing the (Hg^^^^ cis-hexamer [Fig. [5] and 
Fig. IHt] and the corresponding trans-hexamer structure 
[Fig. [6] and Fig. |9]i] where all six C atoms from the 
hexagon decorated with H adatoms share favorable 
tetrahedral coordination. For bigger trans-clusters with 
fully hydrogenated carbon rings [Hio^, His^, Hie^ and 
H22^ in Fig. [7] the H binding is further enhanced, slowly 
approaching the value encountered in graphane. 
A similar investigation has not been carried out for cis- 
clusters since for this class of H adsorption configurations 
trends are less obvious. A full monolayer of H atoms 
can not be adsorbed on the same side of the graphene 
sheet. The maximum H coverage on graphene resulting 
in one or several stable adsorption configurations is also 
unknown. Thus, the trends in the structure and stability 
of bigger H cis-clusters on graphene is an open issue out 
of the scope of the present study. 



A. Trends in stability of large trans-clusters 

To rationalize the DFT results obtained for the big 
trans-clusters. Fig. [71 we construct in the following a sim- 
ple model that reproduces the trends in stability of these 
clusters. In the model we identify closed hydrogenated 
carbon hexagons and a maximal ratio of inside atoms to 
edge atoms in the cluster as the most important struc- 
tural motifs causing high cluster stability. On the other 
hand, removing a H atom from a closed hydrogenated 
carbon hexagon, i.e. introducing an embedded M-dimer, 
is modeled with an energy cost. Graphane, which can 
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O-cis 0-trans H8-cis H8-trans 

Z H-Co-Co 104.84 106.48 Z H-Co-Co 100.81 107.43 

Z H-Co-Ci 101.75 106.18 Z H-Co-Ci 98.31 105.54 

Z Ci-Co-Co 116.01 111.42 Z Co-Co-Co 120.00 108.47 

Z Ci-Co-Ci 113.68 114.56 Z Ci-Co-Co 115.51 113.77 

TABLE V. The angles (in deg.) calculated for dimers (O-cis and 0-trans) and hexamers (He^'^'-cis and (He^^^^-trans) in Fig. 

m 



be considered as an infinite trans-cluster with no edge 
atoms, thus represents the highest achievable cluster sta- 
bility. 

In the model, the binding energy per H, E{model), is 
defined as: 



E{model) = (iV, 



inside-^inside 



NMEM)/n. 

(2) 

where Ninside and N^dge are the number of inside and 
edge atoms, respectively and where Nm is the number 
of carbon sites where the lack of an H atom causes the 
appearance of an embedded M-dimer. The three energy 
terms, Ej^^^^^, ^edge^ ^^"^ are fixed from calculated 
DFT values, i.e. they are not considered ajustable param- 
eters. For fi'j^si^e ^6 use 2.45 eV, which is the H-binding 
energy of graphane when the lattice constant of graphene, 
1.42 A is used (using the self-consistent graphane lattice 
constant of 1.46 A it would be 2.49 eV). For E^^^^ we use 
2.08 eV, which is the H binding energy of 7?g'^^^— trans 
that has H atoms exclusively in edge sites. Finally, the 
energy penalty of introducing M-dimers, Em, we extract 
from the calculated DFT energy difference between the 
Eq^^^^ -trans and the E^s'^^^-trans clusters [see Fig. [6j, 
i.e. Em= 5*2.08 - 5*1.84 = 1.2 eV. 
For each of the clusters from Fig. [7] the relevant param- 
eters and calculated energies are given in Table IVII and 
Fig. [TOl An excellent correlation between DFT val- 



Config. 


^inside 


Nedge 


Nm 


E{model) 
(eV/H) 


E(DFT) 
(eV/H) 




1 


8 


1 


1.99 


2.02 


Hio' 


2 


8 





2.15 


2.15 


Hio'^ 


4 


6 


3 


1.87 


1.83 


H12' 


3 


9 


1 


2.07 


2.11 




4 


8 


1 


2.10 


2.08 


H13' 


4 


9 





2.19 


2.17 


His' 


6 


10 





2.22 


2.22 


H22' 


10 


12 





2.25 


2.23 



TABLE VL The parameters and H binding energies from 
model, E{model), and from the the full DFT calculations, 
E{DFT). 

ues for H binding in trans-clusters from Fig. [7] and the 
results obtained using Eq. [2] demonstrates that the gen- 
eral trends in the stability of graphane-like patches on 
graphene can be rationalized from the simple model with 
the parameters all determined by DFT calculations. 



2.2 

I 2.1 



1.9 



1.8 



0/ 



1.8 1.9 2 2.1 2.2 
E{DFT) (eV/H) 

FIG. 10. The correlation between hydrogen binding energies 
for the clusters in Fig. [71 obtained from DFT calculations 
and those produced applying Eq. (2] 



B. Trends in H diffusion barriers in cis-clusters 

Even in the cis-structure with the highest H binding 
energy, i.e. the configuration Hq, the H binding energy is 
smaller than in the gas-phase H2 molecule. Hence, all cis- 
clusters considered in the present study are metastable 
against associative H2 desorption. Yet, at sufficiently low 
temperatures this process will be hindered due to the 
considerable activation energies that we calculate. Ac- 
cording to our previous studies^SiS^ the lowest activation 
energy for H2 recombination from H dimers is calculated 
for the P-dimer state. In other dimcr structures and for 
many of the large H clusters, the recombination includes 
one or several diffusion steps of H atoms leading to the 
formation of the embedded P-dimer, followed by the II2 
associative desorption. 

To lay the grounds for understanding these events, we 
performed a comprehensive set of calculations regarding 
energy barriers for H diffusion on graphene, consider- 
ing most of the structures presented in Figs. [3][5] Ac- 
cording to the universality concept of N0rskov and co- 
worker o"'^'^" , the activation energies for diffusion of atoms 
on (metal) surfaces correlate with the differences be- 
tween their binding energies in the initial and final states. 
Correlation diagrams-so-called Br0nsted-Evans-Polanyi 
(BEP) plots-reveal that weaker bound adsorbates gener- 
ally diffuse with lower activation energies. The BEP plot 
in Fig. [11] clearly demonstrates this trend for H diffusion 
at graphene. The barrier for diffusion of an isolated H 
atom (AE=0) is 1.14 eV. However, the barrier vanishes if 
the final state is by ^2 eV more favorable than the initial 
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FIG. 11. Br0nsted-Evans-Polanyi plot for selected H cis- 
clusters from Figures 1 2151 For each of the considered H clus- 
ters we calculated the energy barrier (Ea) for H diffusion to 
one of nearest non-hydrogenated C sites. AE is a difference 
in total energies of a cluster configuration produced upon H 
diffusion and the initial one. The dashed lines represent cor- 
responding linear relations between Ea and AE for structures 
with odd and even number of H atoms given by Eq. [S] 



one, or increases above 2 eV for diffusion in the opposite 
direction. For most of the configurations investigated in 
the present study the activation energies for diffusion of 
atomic H are in the range from 0.7 to 1.5 eV. Applying 
the linear least square fitting to the DFT values calcu- 
lated for Eq (the activation energy for H atom diffusion 
between adjacent C sites) and Ai5 (the energy difference 
between configurations after and prior to H diffusion), we 
arrive at the folowing expressions for clusters with n H 
atoms: 

. . J 1.0eV + 0.5A£;, n even . 
J^a[n) - j 1 2eV + 0.5A£;, n odd. ^"^^ 

We did not find clear arguments for observed energy shift 
of 0.2 eV in the between configurations with even 
(dimers, tetramers and hexamers) and odd (trimers and 
pentamers) number of H atoms. Yet, Eq. [3] and the 
BEP plot in Fig. [TT] demonstrate that for H structures 



on graphene a fairly accurate estimate of their stability 
against H diffusion can be obtained from total energies 
of initial and final states, without explicit calculation of 
the corresponding activation energy Ea. 



V. CONCLUSIONS 

We investigated an extensive set of adsorption configu- 
rations of H atom clusters on graphene, the H atoms be- 
ing either all on one side of the graphene (cis-clusters) or 
on both sides in an alternating manner (trans-clusters). 
The binding energy per H atom in general increases with 
the size of the clusters. The value of 0.77, calculated 
for monomers, increases to 1.6 eV in cis- hexamers and 
to 2.08 eV in trans- hexamers. H-H interactions appear 
to favor cis-cluster shapes having H atoms in O- and P- 
positions with respect to each other and to favor trans- 
clusters having H atoms in 0-positions with respect to 
each other. Very stable trans-clusters with 13-22 H atoms 
were identified by optimizing the number of H atoms in 
ortho-trans-positions and thereby the number of closed, 
H-covered carbon hexagons. For such clusters, H binding 
energies up to 2.23 eV were found. For the cis-clusters, 
associative H2 desorption was investigated. It was found 
that desorption occurred from P-dimers embedded in the 
clusters or from pairs of H atoms that first rearranged via 
diffusion into embedded P-dimers, with barriers ranging 
from ^1.4 to '^1.7 eV. Upon formation of clusters with 
an odd number of adsorbates a non-zero spin-density is 
induced. However, for the most stable cis-configurations 
the spin-polarization is rather weak. The most stable cis- 
configurations with an even number of H atoms are all 
non-magnetic. Yet, some of the very stable trans-clusters 
are magnetic, which opens routes for design of magnetic 
materials based on graphene functionalized with hydro- 
gen. 
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